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Local rheological measurements in the granular flow around an intruder 
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The rheological properties of granular matter within a two-dimensional flow around a moving 
disk is investigated experimentally. Using a combination of photoelastic and standard tessellation 
techniques, the strain and stress tensors are estimated at the grain scale in the time-averaged flow 
field around a large disk pulled at constant velocity in an assembly of smaller disks. On the one 
hand, one observes inhomogeneous shear rate and strongly localized shear stress and pressure fields. 

On the other hand, a significant dilation rate, which has the same magnitude as the shear strain 
rate, is reported. Significant deviations are observed with local rheology that justify the need of 
searching for a non-local rheology. 

PACS numbers: 45.70.-n 83.80.Fg 


I. INTRODUCTION 

The description of the flow properties of granular mate¬ 
rial remains a true practical and fundamental challenge. 
This challenge comes from the inherent non linear inter¬ 
actions betweens the grains and from the small separa¬ 
tion of length scales. As a matter of fact, no (attractive) 
force exists between two grains when the contact is bro¬ 
ken and the flow scale is often not much larger than the 
elementary grain scale. Many tools have been developed 
over the past decades to tackle this issue and extract the 
kinematic properties of grain flows, underlying the open 
problem of defining a single rheology for granular matter. 
In the past decade, a local rheology emerges based on a 
unique relation of the friction coefficient /a with the so- 
called inertial number I corresponding to the shear rate 
made dimensionless with confining pressure and density 
p], [H • This rheology works quite well by rescaling various 
experimental and numerical data into a consistent pic¬ 
ture for parallel flows - such as Couette flows or flows 
down inclines - or weakly non parallel flows in a wide 
range of flow rates |1. 3, 4]. Nevertheless, this local des¬ 
cription falls short of describing non-parallel flows where 
the streamlines are far from parallel and also quasi-static 
regimes close to the “liquid-solid” transition 0. Non¬ 
local effects, where the stress not only depends on the 
strain rate but also on its spatial variations, have been re¬ 
cently reported both experimentally and numerically (gi¬ 
ll]. These effects may arise from plastic rearrangements 
inducing long range correlations and cooperative beha¬ 
viors. Non-local modeling for granular flows is thus now 
considered to extend the u(I) rheology by the introduc¬ 
tion of a diffusive term |8l4lll|. This captures well the 
average flow fields reported in some numerical simula¬ 
tions and experiments where the flow is observed to be 
strongly inhomogeneous. Such developments were inspi¬ 
red by kinetic theory, where non-locality arises from the 
diffusion of a granular temperature into a granular flow 

0G1. 

In this paper, we present quantitative measurements 
in the two-dimensional flow of small photoelastic disks 


around a larger intruder disk, where the local strain rate 
and the local stresses are recorded simultaneously. This 
allows one to test the local rheology framework in the 
case of a strongly non-parallel flow. The four invariants of 
the local strain rate tensor and local stress tensor reveal 
that the flow is strongly localized around the intruder. 
A tentative rheological analysis sheds light on a signifi¬ 
cant amount of dilation and some spatial variations of 
the friction coefficient that bring out new challenges for 
modeling. 


II. EXPERIMENTAL SETUP 

The experiments consist in pulling an intruder disk 
at constant velocity Vo into a two-dimensional assembly 
of about 8 x 10 3 smaller photoelastic disks in a set-up 
sketched in Fig. QJa) which is adapted from 0. This 
granular assembly is a bidisperse mixture of disks of dia¬ 
meters 4 and 5 mm to avoid any possible crystallization, 



Figure 1 . (Color online) (a) Sketch of the experimental se¬ 
tup. (b) Snapshot of the intruder moving from left to right 
through the granular packing, (c) Same image obtained bet¬ 
ween crossed polarizers and displaying transient chains of 
contact forces. 
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and made of polyurethane (density p s = 1280 kg/m 3 ) 
which is photoelastic. These disks are placed in bet¬ 
ween two horizontal glass plates with a gap thickness 
of 4 mm, slightly larger than the disk height h = 3.2 mm, 
and closed by lateral walls that form a square area of 
L x xL y = 400x400 mm 2 ~ 90x90 d 2 , where d g = 4.5 mm 
is the mean diameter of the granular assembly. By va¬ 
rying the number of disks, the global packing fraction 
<po of the model granular medium has been varied in the 
range 0.76 < </>o < 0.80 that corresponds to a dense pa¬ 
cking but remains always significantly smaller than the 
jamming point </>j ~ 0.83 [3]. The larger intruder disk 
of diameter d = 16 mm corresponding to about 4 d g is 
pulled in the ^-direction with a steel wire attached to 
a linear stepper motor at a constant velocity Vo in the 
range 10 -4 < Vo < 3 mm s -1 . In order to prevent the 
intruder from tilting and to ensure minimal interaction 
between the wire and the grains, the steel wire is welded 
on the top face of intruder and we designed the intruder 
slightly taller than the disks. As a result, the overhead 
wire hardly sags and has minimal contact with the grain 
assembly. Since the wire is very stiff, the intruder does not 
tilt more than a few degrees. The steel wire is essentially 
inextensible and ensures that the pulling device does not 
deform when the intruder travels across the medium. A 
piezoelectric sensor records the time evolution of the drag 
force F(t) exerted by the granular medium onto the in¬ 
truder disk during its displacement. Once the intruder 
placed in the granular medium at one side of the cell at its 
initial position (Xi,Yi) ~ (—0.4L X ,0) and before its pul¬ 
ling motion to the other side of the cell towards its final 
position ( Xf,Yf ) ~ (0AL x , 0) close to the front wall, the 
whole packing is first homogenized by horizontal vibra¬ 
tions of the bottom plate in the transverse ^/-direction. All 
experiments presented here were performed at constant 
volume with a fixed cell size (constant <po). During the in¬ 
truder displacement, snapshots are recorded using a high 
resolution camera (2048 x 2048 pixels) at a frame rate 
/ such as Vo// < d g / 10. To avoid any possible boun¬ 
dary effect in the results, images and subsequent data 
analysis are restricted to a central square region of in¬ 
terest of size 160 x 160 mm 2 ~ L x j 2.5 x L y / 2.5, leading 
thus to a spatial resolution of about 0.1mm ~ d g / 45. 
Each grain is thus composed by about 1.6 x 10 3 pixels 
in the images. Since the forces at stake are of the order 
of 0.1 N, the in-plane compression is smaller than 0.01%, 
so that we consider a constant diameter in the detection 
algorithm. The granular packing is backlit with a large 
uniform source of linear polarized light and a switchable 
circular analyzer is mounted on a stepper motor synchro¬ 
nized with the picture acquisition. Such a device gives one 
alternated cross-polarized and direct light pictures at the 
effective frame rate of //2. Assuming little change in the 
packing between two successive pictures, one combines 
structural information extracted from the direct light 
pictures [Fig. W[b)] and stresses measurements extracted 
from the cross-polarized pictures [Fig. |T](c)] to quantify 
stress and strain tensors at the grain scale. Considering 


the flow geometry, the results will be presented in polar 
coordinates (r, 6) centered on the intruder, with 0 = 0 in 
the x-direction of the intruder motion. 


III. IMAGE ANALYSIS 


For each direct image, we extract the position of each 
grain including the intruder. By tracking each grain bet¬ 
ween two successive images, we compute its velocity from 
its displacement. To limit noise, the velocity gradient ten¬ 
sor Vv at the grain location is not calculated by differen¬ 
tiation but by integration along a closed contour around 
the grain |Jj5. UK] : 


W=- 


' ® n dZ, 


(i) 


where S is the surface delimited by the correspon¬ 
ding contour, dl is the elementary displacement on this 
contour for which the outward normal unit vector is n, 
and (8) denotes the tensor product. Here, as each grain is 
surrounded by a finite and discrete number of neighbors, 
we compute the integral of Eq. (pQ) along the edges of its 
Voronoi cell. The strain rate tensor G is then deduced 
for each grain by 


G = i(Vv + Vv r ), (2) 

where the exponent T refers to the transpose matrix. The 
strain rate tensor is finally interpolated along a cartesian 
grid of step dx = dy = 0.4 mm ~ d g / 10. 

Another important part of the image processing is the 
computation of the stress tensor cr. From each direct 
image, the Delaunay triangulation and the Voronoi tes¬ 
sellation are constructed to extract the neighbor network 
and contact network of each grain at each time t. From 
the cross polarized image at the same given time t (more 
exactly that only differs by 1 //), we determine then the 
contact forces between grains by calculating the spatial 
gradient of the intensity inside triangular sectors defined 
by one particle center and the corresponding side of the 
Voronoi cell mil. The stress tensor for each grain is then 
computed as fl5l. 1 181 ] 


<T = ^E fi i l8ri J’ ( 3 ) 

1 

where f \j is the contact force between the two grains i and 
j in contact and r ij is the vector of modulus di/2 pointing 
from the center of grain i of diameter di towards the 
center of grain j. The stress tensor a is then interpolated 
on the same Cartesian grid as for the strain rate tensor 
G. Note that the present 2D stress is a force per unit 
length and not a force per unit area as the real stress 
cr/h would be. 

Finally both strain rate and stress tensors G and cr are 
decomposed into an isotropic part and a deviatoric part 
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0 : 


G — el + Gd, 

( 4 ) 

<7 = -pi + er d , 

( 5 ) 


where I is the unit tensor, e = ^tr(G) = |divv is the 
dilation rate of the 2D grain flow, Gd is the shear rate 
tensor, p is the pressure and a& is the shear stress tensor. 
In the following, we shall present the results using only 
scalar quantities corresponding to the invariants of the 
strain rate tensor G and the stress tensor a. The first in¬ 
variants are the dilation rate e and the pressure p and the 
second invariants corresponding to the deviatoric parts of 
the two tensors are the shear rate 7 and the shear stress 
r as named usually in rheology 0 and defined as : 

\ 1/2 

tr(G^)J , (6) 

/I \ 1/2 

r = (-tr(<r’)J . ( 7 ) 



The drag force F on the intruder can be inferred from 
the image analysis described above by the integration of 
the local stress tensor around the intruder : 

F= (f an • e x dl. ( 8 ) 


For the image analysis, we leave out the images where 
the intruder is too close to the edges of the image and 
only focus on the images where the intruder is in a central 
zone (—60 ^ X ^ 60mm, i.e. —L x /1 < X < L x /7). This 
allows one to have local information in the granular me¬ 
dium up to at least 22 mm ~ 5 d g away from the intruder 
center. This should be a large enough zone of analysis as 
the radial extension of the velocity field perturbation is 
expected to be here of about 2 d g + d /4 ~ 13 mm ~ 3 d g 

m 

The time evolution of the drag force F(t) extracted 
from the image analysis [Eq. (0 )] is reported in Fig. [21(a) 
together with the global measurements from the piezoe¬ 
lectric sensor. Both show strong fluctuations that are un¬ 
doubtedly associated with the creation/breaking of force 
chains inside the granular material but no significant slow 
variations. The fluctuations of these two signals are ra¬ 
ther well time-correlated even if their amplitudes are dif¬ 
ferent which may be related to the very different mea¬ 
surement technics. As a matter of fact, the force sensor 
measures not only the drag on the intruder but also the 
friction of the intruder with the bottom and the pos¬ 
sible friction of the wire attached to the intruder. By 
contrast, the photoelastic technique only measures inter¬ 
grain contact forces. Anyway, the absence of slow time 
variations suggests that a stationary regime is reached 
where the force is constant on average, in agreement with 




Figure 2. (Color online) (a) Instantaneous drag force F(t) 

on the intruder obtained either (-) by the force sensor or 

(—) by the image analysis from the photoelastic images for a 
intruder of diameter d — 16 mm in motion at the velocity Vo = 
10“ 2 mm s -1 in a 2D granular medium of packing fraction 
0o = 0.76. (b) Time-averaged drag force (F) as a function of 
Vo. Data from (v) force sensor and (a) photoelastic images 
together with (—) the mean value (F) = 0.22 N averaged over 
the explored velocity range 10 -5 < Vo < 10 mm/s. 


previous observations for the motion of objects in a gra¬ 
nular material 01 SHU- In the following, we only focus 
on time-averaged quantities along the intruder motion 
in the central part of the cell. Figure EJb) shows that 
the mean drag forces (F) obtained from the force sensor 
and from the photoelastic signal are very close to each 
other whatever the pulling velocity Vo, which gives confi¬ 
dence in our local image analysis. For Vo ranging within 
about five decades, ( F(t)) is roughly constant, in agree¬ 
ment with the experiments mentioned above. 

In the next section, we will focus on the time-averaged 
local quantities that are extracted from image analysis of 
both direct and polarized images. These different quan¬ 
tities characterizing the granular flow will be normalized 
by the natural length scale d g and time scale d g /Vo. As 
there is no external stress applied on the lateral boun¬ 
daries and no “hydrostatic pressure” from gravity in our 
configuration, the typical stiffness k = (n/A)Eh of the 
contact between two cylinders jUj] will be taken as the 
stress scale, where E is the effective Young modulus of 
the disks of height h with here the value k ~ IN mm -1 . 


IV. LOCAL MEASUREMENTS 

A typical velocity field of the granular medium around 
the moving intruder is displayed in Fig.[3](a). As expected 
from previous kinematic description pjj [ 21 ] , the velocity 
perturbation induced by the intruder motion is found to 
be strongly localized around it. The corresponding local 
packing fraction 0 is displayed in Fig.[3](b) : 0 is observed 
to be quite homogeneous in space and close to the global 
value 0 o except in a very narrow crown upstream of the 
intruder and in an elongated triangular wake downstream 
where 0 is significantly smaller. This downstream wake 
has been reported in detail by jib] in similar 2D experi¬ 
ments and follows directly from the absence of confining 
pressure on the cell boundaries or gravity, as the present 
setup is horizontal. 
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Figure 3. (Color online) (a) Velocity field around an intruder 
disk of diameter d — 16 mm moving from left to right at the 
velocity Vo = 2.7 mm s -1 into a granular medium of global 
packing fraction 0o = 0.76. Vectors indicate the grain velo¬ 
city direction and the colormap encodes the modulus of the 
dimensionless velocity v/Vo. (b) Corresponding local packing 
fraction 0. 
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Figure 4. (a) Local packing fraction 0 as a function of the 
radial position r/d 9 from the intruder along the azimuthal 
direction 0 = (10 d= 3)°. (b) Local packing fraction 0 as a 
function of the azimuthal position 6 for the radial position 
r — (11 =b 1) mm (i.e. r/d 9 ~ 2.4 zb 0.2 from the intruder). 
Data are obtained for d = 16 mm, Vo = 10 -2 mm s _1 , and 
0o = 0.76. 


The radial and azimuthal profiles of 0 are displayed 
in Figs. Ufa) and (b) along one azimuthal or radial po¬ 
sition respectively. The data points have been obtained 
by digitizing the 2D interpolated grids onto annular sec¬ 
tors (r — dr, r + dr ), (0 — dO, 0 + d0), with dr = d g /S and 
dO = 7t/ 32 Figure Ufa) shows that 0 is quasi-constant 
and very close to 0o = 0.76 upstream of the intruder 
except in a very narrow zone of size one grain diameter 
close to the intruder where 0 decreases significantly and 
abruptly with r. The azimuthal profile close to the intru¬ 
der [Fig. 01(b)] shows that there is no significant varia¬ 
tion of 0 around the intruder in a large azimuthal range 
(—7r/2 < 6 < 7r/2) except in the downstream wake ob¬ 
served in Fig. [3Kb). 

The maps of the four scalar invariants i, 7 , p and r 
extracted from the strain rate tensor and stress tensor 
respectively are displayed in Fig. O All these quantities 
vary significantly upstream of the intruder with some dif¬ 
ferences. Concerning the strain rate field [Fig. EJa)], one 
observes a zone of high positive s extending about one 
intruder diameter upstream and in a quite large azimu¬ 
thal extension (— tt/2 < 0 < 7 r/ 2 ), with thus a “banana 
shape”. The shear rate 7 in Fig. [5Kb) displays a similar 



Figure 5. (Color online) Maps of the four dimensionless sca¬ 
lar invariants of the (a,b) strain rate tensor and (c,d) stress 
tensor around an intruder of diameter d — 16 mm moving 
from left to right at the velocity Vo = 2.7 mm s -1 through a 
granular packing of global packing fraction 0o = 0.76 : (a) di¬ 
lation rate ed 9 /Vo , (b) shear rate jdg/Vo , (c) pressure p//c, 
and (d) shear stress r/k. 


banana shape with a radial upstream extension of about 
half the intruder diameter and a similar azimuthal exten¬ 
sion. This correlation between e and 7 is related to the 
fact that dense granular media are known to dilate upon 
transient shear |27l. 1281] . Finally the stress fields [Figs.EKc- 
d)] differ from the strain fields. A zone of high pressure 
p elongated in the upstream direction appears in front of 
the intruder and a zone of high shear stress r of similar 
shape is also observed. There is thus a clear link between 
p and r. The relation between stress and strain appears 
however less clear from these four maps. 

To quantify more precisely the spatial variations of the 
strain and stress tensors, the radial and azimuthal pro¬ 
files of the normalized four scalar invariants are displayed 
in Fig.[ 6 ]for the same global packing fraction (0o = 0.76) 
but three intruder velocities Vo ranging over 4 decades. 
As previously, the data points have been obtained by 
digitizing the 2D interpolated grids onto annular sec¬ 
tors (r — dr, r + dr ), (0 — dO, 0 + d0), with dr = d g /8 
and dO = 7r/32. The rather good collapse observed in 
Figs. EKa-d) of the three curves for the three veloci¬ 
ties shows that i and 7 scales with Vo as expected. By 
contrast, Figs. EKe-g) confirm that p and r do not de¬ 
pend significantly on Vo which is consistent with the fact 
that the measured drag force (F) does not depend si¬ 
gnificantly on Vo in Fig. [2Kb). In the lin-log plots of 
Figs. EKa,c,e,g), the strong monotonic decrease of i, 7 
and p with increasing radial distance from the intruder 
is in agreement with the exponential radial decrease al¬ 
ready reported for 7 |]]|. By contrast, r displays a non- 
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Figure 6. Variations with the normalized radial distance r/d g 
along the azimuthal direction 6 m (10 d= 3)°, and with 6 
for r — (11 ± 1) mm of the four dimensionless scalar inva¬ 
riants : (a,b) dilation rate ed g /V o, (c,d) shear rate jd g /Vo, 
(e,f) pressure p//c, and (g,h) shear stress r/k , for 0o = 0.76 
and Vo = 10 -4 (x), 10 -2 (o), and 1 mm s -1 (+). (—) Cosine 
function as a guideline for the eye. The grey zone corresponds 
to the intruder extension up to r = d/2 — 8 mm. 


monotonic evolution with a first increasing and then de¬ 
creasing evolution with radial increasing distance from 
the intruder surface. Figures [ 6 jb,d,f,h) show that the four 
invariants have a maximal value in the direction of mo¬ 
tion (0 = 0) but quite different variation with 6. The 
variations of e and 7 are close to a cosine shape. Such 
a cosine shape has been already observed for 7 in the 
granular flow around a cylinder in vertical penetration 
motions 0,113 which is thus a consistent result. The 0 
variations of p and r are in contrast far from a cosine. 
The present observed strong localization of p and r with 
both r and 0 may explain the very small length scale re¬ 
ported previously for the force interaction of a moving 
object with both bottom (H] or lateral [29j walls. Note 
that for the lower speed experiment (Vo = 10 - 4 mm/s) 
the radial variations of e and r in Figs.[ 6 ](a) and (g) are 
quite different from the others, which will be discussed 
in the next Section. 




Figure 7. (a) Dimensionless dilation rate ed g /\ b as a function 
of shear rate 'ydg/Vo and (b) dimensionless shear stress r/k 
as a function of pressure p/k for an intruder of diameter d — 
16 mm. Same data points as in Figs.[6ja, c, e, g) together with 
(A) for 0o — 0.80 and Vo = 10 -2 mm s -1 . 


V. RHEOLOGICAL CURVES 

Let us now look at the relations between the four sca¬ 
lar invariants £, 7 , p, and r to search for the existence of 
a possible local constitutive law for the granular medium 
in the present non-parallel flow around a moving intru¬ 
der. Despite some dispersion, Fig. [T^sl) shows that the 
dilation rate e increases with the shear rate 7 . This cor¬ 
relation between e and 7 is reminiscent of the so-called 
Reynolds dilatancy [30], where dilation usually increases 
quadratically with shear strain 6 ~ 7 2 ( 27 ). With also 
some dispersion, Fig. [7Kb) shows that the shear stress r 
increases with the pressure p which justifies the intro¬ 
duction and use of a friction coefficient. Note that we 
ensure that the principal directions of the shear rate and 
shear stress tensors are nearly parallel in the upstream 
region where a local dilatancy coefficient D and a fric¬ 
tion coefficient p will be computed, in agreement with 
the requirements pointed out by [!j. 

In order to test the possible local rheology, we now look 
at both relations of the dilatancy coefficient D = e /7 and 
the friction coefficient p = r/p with the inertial number 
I = 7 d g \Jps4>/p which is the ratio of the inertial time 
scale d g /^/p/p s </> at the grain scale and the flow time 
scale I /7 , and which has been shown to be the rele¬ 
vant dimensionless number in simple parallel flows GM3- 
The data shown in Figs. [ 8 ] and [9] are the D, p and I va¬ 
lues averaged over the radial range 2 ^ r/d g < 6 along 
different azimuthal directions, from the motion direction 
near 0 = 0 [Figs.[ 8 Ka) and [51(a)] towards the equatorial di¬ 
rection near 0 = tt/2 [Figs.[ 8 ](d) and [51(d)]. One observes 
that D does not depend significantly on / with about the 
value D = 2 db 1 that remains about the same whatever 
the azimuthal direction 0. Concerning the local friction 
coefficient, one does not observe any significant increase 
of p with /, which is due to the very low / values explored 
here corresponding to a quasi-static regime (/ < 10 -2 ). 
In some directions [<9 < 45 0 in Figs. [51(a) and (b)], one 
may even observe a significant decrease of p with increa¬ 
sing I at very small values (/ < 10 -5 ). This is linked to 
the differences reported in Fig. [ 6 ] for the lowest speed ex¬ 
periments (lowest I). This behavior is expected to arise 
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Figure 8. Local dilatancy coefficient D as a function of the 
inertial number I with an average over the radial extension 
r/dg < 6 along four azimuthal directions : (a) 6 = 10 d= 3°, 
(b) 6> = 45±3°, (c) <9 = 60 ±3°, and (d) <9 = 80 ±3°. Same 
data symbols as in Fig. [6] together with Vo — 10 -3 (□), 10 -1 
(o), and 2.7mm s -1 (v), for 0 O = 0.76. 



Figure 9. Local friction coefficient /i as a function of the 
inertial number I with an average over the radial extension 
r/dg < 6 along four azimuthal directions : (a) 6 — 10 ± 3°, 
(b) 0 = 45 ± 3 °, (c) 0 = 60 ± 3 °, and (d) 0 = 80 ± 3 °. Same 
data symbols as in Fig. [8] 

from theoretical models m but was not yet reported ex¬ 
perimentally, and may explain the stick-slip like behavior 
observed for intruder dragged at very low velocities (23| . 

The key point is however the significant dependence 
observed for // with the azimuthal direction 0. Indeed, a 
significant decrease by a factor about 3 is observed from 
the direction of motion where (i ~ 0.15 [Fig. [9^a) for 
6 ~ 10 °] towards the equatorial direction where (i ~ 0.05 
[Fig. (9](d) for 0 ~ 80°]. This strong 9 variation of //(/) 
means that such a local rheology fails here to describe 
the entire granular flow around an intruder. This experi¬ 


mental result is consistent with recent theoretical findings 
0 where the friction coefficient is found larger in front 
and behind the intruder and smaller on the sides. Our 
present result confirms experimentally these theoretical 
predictions and strongly suggests that non-local rheology 
is thus needed to fully describe non-parallel flows. 


VI. CONCLUSION 


The present paper demonstrates the possibility of com¬ 
putation of the local stress tensor together with the lo¬ 
cal strain tensor in a two-dimensional granular packing 
of photoelastic disks from a detailed analysis of cross- 
polarized images. This possibility is here used to test 
a possible local rheology in a strongly non-parallel flow 
consisting in the grain flow around a larger intruder disk. 
The analysis of the spatial variations of the four inva¬ 
riants of the strain rate tensor and stress tensor reveals 
that the flow is strongly localized around the intruder and 
that the stress is very localized in the direction of mo¬ 
tion. A positive dilation rate increasing roughly linearly 
with the shear rate is put in light and allows one for the 
determination of a dilatancy coefficient D for the present 
stationary non-parallel granular flow. This dilatancy co¬ 
efficient is essentially independent of the azimuthal direc¬ 
tion 0. By contrast, the friction coefficient /i that relates 
shear to normal stresses shows significant 6 dependence. 
This spatial dependence demonstrates that local rheo¬ 
logy may be not sufficient to describe strong non-parallel 
flows such as the present flow around a cylinder. While 
such spatial variations of fi have been recentl y p redicted 
by non-local modeling in intruder geometry [lOj, the fi¬ 
nite level of dilatancy D of stationary flows has not been 
reported yet theoretically or numerically. The present fin¬ 
ding may thus put further constraints on the formula¬ 
tion of continuum models by opening up novel challenge 
in modeling and understanding non-parallel compressible 
granular flows using compressible non-local rheology. Ex¬ 
perimentally, checking for non-local rheology is challen¬ 
ging as the computation of Laplacian terms - at the core 
of these non-local models [8|, Hq( ~ is hampered by finite 
resolution and large fluctuations. 

Similar experiments with an imposed vibration of the 
packing could be developed in order to uncouple the effect 
of the inhomogeneous excitement created by the pulling 
of the intruder and the effect of homogeneous excitement 
( 32 M 35 II . These studies should help to a better understan¬ 
ding of the drag force with depth in the case of the verti¬ 
cal motion under gravity 0, [37[ and the direction of the 
motion in the penetration-extraction problem [38|. [39| in 
order to be extended for the understanding of the animal 
locomotion in sand [40]. 
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